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i-^ \ We consider a quantum particle constrained to a curved layer of a 

" ' constant width built over an infinite smooth surface. We suppose that 

Ch ' the latter is a locally deformed plane and that the layer has the hard- 

en ' 



wall boundary. Under this assumptions we prove that the particle 



Q^. Hamiltonian possesses geometrically induced bound states. 

>'• 

^ : 1 Introduction 

Relations between geometry and spectral properties are one of the trademark 
topic of mathematical physics, in fact, an abstraction of various acoustic 
problems can be found in the roots of this discipline. For a long time, how- 
ever, it seemed that such questions were restricted to the area of classical 
physics. This was mostly because geometrical properties of quantum sys- 
tems were supposed to result from their dynamics, and as such they were not 
accessible to experimenters choice and manipulation. Interesting results de- 
pendent on geometry existed, of course, such as chaotic behaviour of certain 
quantum billiards observed first in ||BGS|] and studied in numerous subse- 



quent papers, but they remained to be mostly academic exercises. 



The situation has changed dramatically with the advent of mesoscopic 
techniques which allow us to produce tiny structures of various shapes de- 
vised and reproducible in the laboratory and yet small enough to exhibit 
quantum effects. Moreover, the physical nature of these objects makes it 
possible to describe them by a simple model in which a free particle (with an 
effective mass) is confined to the spatial region of the structure - see ||Dli|| and 



references therein - borrowing the terminology one can speak about quantum 
waveguides, resonators, etc. 

Theoretical studies of such systems have brought interesting results, some 
of them being purely quantum, without a classical counterpart. Among 
the most beautiful is the binding effect of curvature due to which infinitely 
extended regions with hard walls and a constant width can exhibit local- 
ized states. The effect was in fact observed a long time ago - see, e.g., 
d(Jl| , |d(J2| , [lbl| | and references therein - in formal attempts to justify quan- 



tization on nontrivial manifolds, but only in Pq| it was placed into a proper 



quantum- waveguide context and the existence of geometrically induced dis- 
crete spectrum for a curved planar strip was rigorously proven. This was 
followed by numerous other studies in which the results were improved and 
properties of the bound states were investigated - see, e.g., [pE| , |GJ| , [RB|| . 

Much less is known about analogous system in higher dimensions starting 

from the physically interesting case of a curved layer. This may seem strange 

at a glance, since the leading term of the effective potential for a general n- 

dimensional manifold in R'""^" was computed more than two decades ago 

Toll . However, going beyond the formal limit of infinitely thin layer one has 



to be able to estimate the next terms which is not an easy task. The aim of 
the present paper is to stimulate an investigation of the "two-in-three" case; 
we will concentrate at the simplest case where the deformation of a planar 
layer is infinitely smooth and compactly supported. 

Let us describe the contents of the paper. In the next section we will 
formulate the problem and introduce a technique to handle it based on a 
suitable change of coordinates. The main results are given in Section 3. 
We will show first that under our assumption the essential spectrum starts 
at the first transverse eigenvalue. Then we will present a variational argu- 
ment showing that a local deformation of the layer pushes the bottom of 
the spectrum below this value inducing thus a nonempty discrete spectrum. 
Properties of these bound states will be discussed elsewhere. 



2 Formulation of the problem 

Let So be an open set in M^; its points will be denoted by g = (g^, q^) G M^. 
Let a regular and simple surface S of class C°° in the space M^ be given by 
a mapping 

p:So^M=^:{g^p(g)GE} (2.1) 

such that the vectors p^^ = d^p := dp/dq^, // = 1, 2, are linearly independent. 
We have in mind in this paper surfaces diffeomorphic to the plane, but since 
we will use different parametrizations, it is reasonable to consider Eg generally 
as a subset of R^ . Under the linear independence condition a unit normal of 
the surface 

3.1 ^ ^ ._ P.i X P,2 



n-.Eo^R'-.^q^ n{q) := ,^'^ ^' , G M.' } (2.2) 

is a well-defined smooth function, which defines an orientation of S. Together 
they determine a layer f2 of a width d = 2a > over the surface S by virtue 
of the mapping 

: ^0 -^ K^ : {(g, m) ^ 0(g, u) := p(g) + Mn(g) G fi}, (2.3) 

where f2o := Sq x (—a, a). 

2.1 Properties of the reference surface 

Recall first basic facts about the three fundamental forms of S. The coef- 
ficients of the first fundamental form I of the surface can be identified with 
the covariant components of its metric tensor 

QtJiu ■= P,t^ ■ P,u 9 ■■= det(5f^^) (2.4) 

while for the second one, II, we use the common notation 

hf,u ■= -n^^.p^u h:= det{h^^). (2.5) 

The Weingarten map h'^ {cf. ||Kli| , Def. 3.3.4 & Prop. 3.5.5]) determines the 
Gauss curvature of S, K := det{h'^) = -, and its mean curvature, M : = 
|Tr(/ip = \g^'^hfj_y. The third fundamental form III := {n^^^.n^i,) may be 
expressed by means of the first and second fundamental forms as follows: 



Proposition 2.1 n^^ . n^^ = -Kg^^ + 2Mh^^ = h^p gP" h„^ 

Proof: The first relation is equivalent to the identity III — 2MII + i^I = 
- c/. ||Kli| , Prop. 3.5.6] or ||LR| , Problem 7.53] for a general dimension. To 
prove the second one, we use the fact that K, M are determined by the 
characteristic equation 

k\ - 2Mk± + K = 0, (2.6) 

where fc+, A;_ are the principal curvatures, i.e. the eigenvalues of hj^. Let T^ 
be the principal direction, i.e. the eigenvector of hj^, corresponding to k±. 
Multiplying {^ by this vector we get hP^h'^^T^ - 2MhJ^T^ + Kd^T^ = 0. 
Since the principal directions forms locally an orthogonal basis in M^, the 
same must hold in the matrix sense, Mh", — 2MhZ + K6'i', = 0. The desired 
equality is then obtained by multiplication with g^,y. | 

Remark 2.2 We use the standard summation convention about repeated in- 
dices; the Greek and Latin ones run through 1, 2 and 1, 2, 3, respectively. The 
indices are associated with the above coordinates by (1,2,3) ^^ {q^,q^,u). 
Furthermore, upper and lower index denote components of contravariant and 
covariant tensors, respectively. Indices are raised and lowered by the corre- 
sponding metric tensor. For instance, the matrix of the Weingerten map is 
given by h'^ = h^pg^". The same applies to the metric tensor G^j, in the 
layer Vt which we shall introduce below. 

Next we use the Jacobian (72 to write down the invariant surface element 

dS := g^(fq = g^dq^dq^. 

It makes it possible define a global quantity characterizing S, namely the 
total curvature of our surface as 



Tot(S) := / Kdi:. 

Suppose that ^ C S is a region encircled by a simple closed curve C of 
class C , then the Gauss-Bonnet theorem |[LR| , 7.6.45] claims that 

Tot(^) + i kgdi = 271, (2.7) 

where kg is the geodesic curvature of C (traversed in the positive sense) 
and i denotes its arc length. For the purpose of this paper it is important 
that the geodesic curvature of a circle in the plane is equal to its reciprocal 
radius. Consequently, an infinite surface obtained by a compactly supported 
deformation of a plane has Tot(S) = 0. 



2.2 Metric properties of the layer 



It is clear from the definition ( p.3| ) that the metric tensor of the layer (as a 
manifold with a boundary in M^) is of the following form 

G G 

G,, :=0,,.0,, (G,,)=| Gi2 G22 I , (2.8) 

1 



where G^jy = g^i, — 2u h^^, + u^ n^^ . n^p. In view of Proposition |2.1j we can 
rewrite the last expression as 

G,,=g,,{5P^-uK){5:-uK), (2.9) 

which makes it easy to compute the determinant because the matrix of the 
Weingarten map h^ has the principal curvatures A;+, k^ as eigenvalues. Hence 

G := det(G^,^) = ^ [(1 - uk+){l - uk^)f = g{l - 2Mu + Ku^f , (2.10) 

where in the second step we employed the relations K = k^k- and M = 
|(A;+ + k_). As above, 

dQ := G 2(1^0 = G2(fqdu 

defines the volume element of the layer. 

The "straightening" transformation employed below requires that the 
mapping defining the layer is a diffeomorphism. In view of the regularity 
assumptions imposed on S and the inverse function theorem it is sufficient 
that G^i, has an inverse bounded uniformly in Q. This imposes a restriction 
of the layer thickness d. Define Pm '■= (niax{||A;j,(g)||oo : z^ = ±, g G Sq})^ . 



It follows from (|2.9|) that G^,^ can be estimated by the surface metric, 

C-g^u < G^, < G+g^,, (2.11) 

where the constants G± := (1 ± ap^) are well defined since pm > by 
definition. Hence for a smooth surface S the bijectivity of is ensured as 
long as a < Pm- 



2.3 Various expressions of the Hamiltonian 

After these preliminaries let us define the Hamiltonian H of our model. As 
we have said, the particle is supposed to be free within Vt and the boundary 
of the layer is a hard wall, i.e., the wavefunctions should satisfy the Dirichlet 
boundary condition there. For the sake of simplicity we set the Planck's 
constant h= 1 and the effective mass of the electron m^, = 1/2; then H can 
be identified with the Dirichlet Laplacian 

H:=-A% on L^{n), (2.12) 

which is defined for an open set fi C M"^ as the Friedrichs extension of the 
operator —A with the domain C^i^Q) - cf. [pE| , Sec. XIII. 15]. 

A natural way to investigate the operator ( |2.12| ) is to pass to the intrinsic 



coordinates (g, u) in which it acquires the Laplace-Beltrami form, 

H:=-G~'^diG^G'^dj on L\Qo,G^d\du). (2.13) 



To find its action explicitly, we employ ( p. 81 ) together with the expression of 
the determinant ( p.lO| ). Then H splits into a sum of two parts, 

given by 

H, = -G-'^d^G^G^''d, = ~d^G^''d,-^G-'G,^G^''d, (2.14) 
H2 = -G~'^dsG^^d3 = -d^.-2 ^""-^ Q (2.15) 

The above coordinate change is nothing else than the unitary transformation 

IJ : L^{VL) ^L^{VLQ,G^d^qdu) : {i> ^ U ^ := i) o (j)} 

which relates the two operators hy H = UHU^^ . 

At the same time, it is useful to have an alternative form of the Hamil- 
tonian which is symmetric w.r.t. G and has the Jacobian removed from the 
inner product. It is obtained by another unitary transformation, 

U : L^{VL) -^ L^i^o) : {i/j ^ Uip := G^ o 0}. 



which leads to the unitarily equivalent operator 

H := UHU-^ = -G--^d^G^^G^"'d^G-'^ - G-'^dsG^d^G-'^ (2.16) 

with the domain 

D{H) := {^ e Wl{no)\ Vg G So : ^(g, -a) = ^{q, a) = 0}, , (2.17) 



where 1^2(^0) is the appropriate local Sobolev space in the sense of p5 
Sec. XIII. 14]. Commuting G~^ with the gradient components we cast the 
operator (|2.16| ) into a form which has a simpler kinetic part but contains an 
effective potential, 

H = -d^G^"'d,-dl + V, (2.18) 

with 

V = F} + F,F\ F, :=(lnG^),,. 

This expression of the potential is valid for any smooth metric Gij. If we 
employ the particular form ( |2.8| ) of the metric tensor, we can write again 
(|2.18| ) as a sum of two parts, H = Hi + H2 with \^ = Vi + V2, where 

Hi = -d^G^^d. + Vi (2.19) 

Vi := ^d^G-'G^''G,, + ^G-^G,^G^''G,, 

= -^G-2G,^G^'^G,, + ^G-'G^'^G,^, + ^G-'G,^G^''^, (2.20) 
H2 = -dl + V2 (2.21) 

Vo := (2.22) 

' {l-2Mu + Ku^y ^ ' 

2.4 Coordinate decoupling 

While the operator Hi + V2 depends on all the three coordinates, in thin 
layers its "leading term" depend on the longitudinal coordinates q only. The 
transverse coordinate u = q^ is isolated in H2 — V2 = — <9|, so up to higher- 
order terms in a the Hamiltonian decouples into a sum of the operators 



H, = 


- -g~--d,g-^g^'d,g-^+K-M\ 


(2.23) 


Hu = 


= -dl 


(2.24) 



This observation is behind the formal hniit a -^ 0+ mentioned in the in- 
troduction [ |dCl| , |dC2| , [Tolll in which the transverse part is thrown away and 
the thin-layer Hamiltonian is replaced by the surface operator Hq, with the 
energy appropriately renormalized. This procedure can be given meaning in 



the perturbation-theory framework, in analogy with ||DE|] , as we shall discuss 
elsewhere. 

Here we use it for motivation purposes. The effective "surface potential" 
K — M ^ can be rewritten by means of the principal curvatures of S as follows 

K-M^ = --{k+-k^y. (2.25) 

If the curvature vanishes at large distances, we get a potential well which 
could imply existence of bound states that would persist in layers of a finite 
thickness. In distinction to the "one-in-two" case of a curved planar strip 
the effective potential may vanish if the surface is locally spherical, fc+ = /c_, 
however, this cannot happen everywhere at a locally deformed plane. Let us 
also remark that similar Laplace-Beltrami operators penalized by a quadratic 
function of the curvature lead on compact surfaces to interesting isoperimetric 
problems Pa], [HL| , |EHq . 



In the next section we shall also need the eigenf unctions {Xn}^=i of the 
transverse operator Hu- They are given by 

{-i/|cosfi;„u if n is odd 
(2.26) 
-i/|sinfi;„M if n is even 

and the corresponding eigenvalues are k^ = (kiTi)'^ with ki = ^. 



3 Spectrum of locally deformed planar layers 

In what follows we shall consider a class of layers over surfaces which are 
smooth local deformations of a plane. More specifically, suppose that the 
deformed part of the surface is ^ C S and denote supp K U supp M = 
p^^{A) =: ^o; it is clear that S \ ^ is a plane with a "hole", not necessarily 
a simply connected one. 

Let (X, 6^u) be a natural representation of S\^ in Cartesian coordinates 
{x^jx"^) G X C M^ given by an isometry 

C : Eo \ A ^ X : {g ^ {x\ x^) =: C{q)}. 



In view of the compact-support assumption we can choose ro > in such a 
way that Br^ := {w G S : \Cop^^(^w)\ < Tq} contains ^ and thus T,^^ := S\i3ro 
is the undeformed plane with the disc of radius ro removed. It is useful to 
introduce a polar-coordinate parametrization of S^q given by the isometry 

V : (ro,oo) X S^ ^ X : {(r,^) ^ {x\x^) =: (r cos^,r sin^)}; (3.1) 

the corresponding metric tensor acquires then the form diag(l,r^). 

3.1 The essential spectrum 

In a planar layer the essential spectrum starts from the lowest transverse 
eigenvalue. We will use the standard bracketing argument ||KIJ|, Sec. XIII. 15] 



in combination with the minimax principle to prove that the same remains 
true after a compactly supported deformation. 

Proposition 3.1 crcss(-^) = [f^ijOo). 

Proof: We cut the layer Q perpendicularly at the boundary of Br^ and 
impose there the Neumann or Dirichlet condition respectively; this enables 
us to squeeze H between a pair of operators 

Hint ® -^erf ^ H < H^^^ © -f^esi' (3-2) 

which have both the form of an orthogonal sum. The spectrum of the interior 
parts is purely discrete, so the essential components are determined by the 
exterior part only, (t^^^{H^^^ © i/f^J = a^ssiH^^t) = a(i7f,j), f3 = N,D. The 
latter can be simply localized employing the polar-coordinate parametriza- 



tion ( |3.1| ) of Srg. In particular, let 



Ur : L^ {Qo \ p'\^ro) X i-a, a)) ^ L^ {{ro, oo) x S^ x (-a, a)) 

be the substitution-type unitary operator {U-pip){r,{},u) := ip{V~^ oC{q),u). 
It is clear that the spectrum of the corresponding exterior Hamiltonians 

contains all points ^^ + 6/3, where e/3 belongs to the spectrum of the s-wave 
radial part, /i^ := —df. — (4r^)^^ in 2y^(ro,oo) with the appropriate b.c. at 
r = tq. We have 

-{d'^)^-^<h'' and h''<-{d^.)n 
4rg 



in the sense of quadratic forms and inf o"ess( — ((5^)/3) = 0. Since tq can be cho- 
sen arbitrary large, the claim follows from (|3.2| ) by the minimax principle. | 

3.2 Existence of Bound States 

Now comes the main result of this paper. We are going to show that the 
conjecture about existence of a discrete spectrum in locally curved layers 
formulated in Sec. p.4| is true, even for layers which may not be thin. The 
variational proof of the following results is based on the idea adapted from 
IGJ, see also |Dg, Thm. 2.1]. 



Theorem 3.2 Suppose that the layer is not planar and the deformation sat- 
isfies the smoothness and compact support assumptions. Then inf cr(if) < k\. 

Proof: Denote the norm in L'^{ilQ,G2d'^qdu) as || ■ Hg*; then it follows 
from ( p^.l3D that the quadratic form associated with our Hamiltonian H is 
given by 

q[^/J]■.= \\H-2'l|;\\l = q,[^jJ] + q^[^|J] 

where 

gi[^] := ll^fV^II^ = (V^,^,GiG^>,,) (3.3) 

g2[^] := lli^l^ll^ = ||G^,3|r. (3.4) 

It acts on Q{H), the quadratic form domain of H. In order to prove the 
claim it is sufficient to find a trial function ip G Q{H) such that 

t[ij]:=q[^P]-Kl\ml<0. 

(a) We begin the construction of a trial function with ip{q,u) := ip{q)xi{u), 
where xi is the lowest transverse-mode function ( |2.26| ) and (p is a. function 



from the Schwartz space iS(M ), arbitrary for a moment. It yields 

giM = [v,„{G^2G^^\xi?)u^,u)^ (3.5) 

[ip,{G^x'i?)uv)^ (3.6) 

[^,{G-2\xi?)u^) (3.7) 



H 



2 _ / iry^\ |2\ 
G - 
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where (■)„ means a "transverse" expectation and the subscripts g, u mark 
the fact that we integrate w.r.t. the corresponding coordinate only. 

Taking into account the exphcit expression ( p.lO| ) for G and using the 



trivial fact that |xi|, \Xi\ ^ire even functions and that we integrate over a 
symmetric interval {—a, a), and consequently, that we can consider just the 
even powers of m in (■), we get 



li/iOv 0/^^l| lOv ^^1 



{G-^lxWlu- f^i{G-^\Xinu = Kg-^ 



we have employed at that the identity ("U^dXip— i^i\Xi\'^))u = 1- By virtue 
of ( p.ll| ), we can estimate the remaining term as 

qi['iP]<C+(ip,^,gh^''V>,u) . (3.8) 

Suppose now that (p{q) = 1 on p^^{Bro) and that the function is radially 
symmetric in the sense (p{r, i)) = (p{r), where ip := ipo C^^ o V. Passing then 
to the polar coordinates {g2(Pq = rdrdd) in ( |3.8| ) we arrive at 

QiWl < C+ / \(p\^r drdd =: C+ ||<^||7>, 

where the dot denotes the derivative w.r.t. r. The r.h.s. of this inequality 
depends on the surface geometry through the constant C+ only. Summing 
up the results we have 

m<C+\\(p\\l + {^,Kg"2^)^. (3.9) 

(b) In the next step we shall specify further the function (^ in a way which 
allows us to make the r.h.s. of ( |3.9| ) arbitrary small. Let us define the family 
{(pa '■ o G (0,1]} by an external scaling (in the region r > tq) of a suit- 
able function. The idea is analogous to |pJ|| or |PE| , Thm. 2.1], however. 



since we deal with a two-dimensional integral we have to be more careful 
about the decay properties. We can adopt for this purpose the mollifier em- 
ployed in ||EV| , PCEZ]| , which is expressed in terms of Macdonald functions 
(or modified Bessel functions in the terminology of | |AS| ]) as 

ip„[r) := mm <^ 1, 



^o(o-ro) 
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Since Kq is strictly decreasing, the corresponding ipa '■= faXi "will not be 
smooth at r = To but it remains continuous, hence it is an admissible trial 
function as an element of Q{H). 

To estimate the first term at the r.h.s. of ( p.9| ), let us compute the norm 
of the scaled function using [|AS| , Sec. 9.6] and [|GR| , 5.54]: 



WaWv 



27r 






Ko(arY r dr 



2n 



^oKo)2 



Ki{tftdt 



aro 






aro 



In aro 



— {crroY In aro — 2aro In aro 



+ (aro) In aro 



K,{aro) 
Ko{aro) 

^iKo)Y 



Ko{aro)j 

Next we use the small- argument asymptotic expressions [|AS| , Sec. 9.6] 

Ko{x 
Ki(x 



-lna; + C>(l) 
+ 0(lnx) 



X 



Ki(x) 



which imply that a; In x ^^ remains bounded as x 



IVaWv < 



0+, hence 



(3.10) 



In aro | 

holds for a positive constant b and aro small enough. 

(c) To handle the second term at the r.h.s. of ( |3.9| ) for ip = ip„ we employ 
the dominated convergence theorem: since \(pcr\ < 1 and y^o- ^ 1 pointwise 
as a — i> 0+, we have 



i^a,Kg2ip^) 



'l,Kg-^),^TotiE) 



by definition. Notice that the total curvature integral is well defined because 
the Gauss curvature K of a surface obtained by a smooth compactly sup- 
ported deformation of a plane belongs to L^{T.o,g2d'^q). In view of 
( |3.1CI| ) and the last formula we have therefore 



t[i;^] ^ Tot(S) 



(3.11) 



12 



as (T -^ 0+. Recall that in Sec. \2,.1\ we have used the Gauss-Bonnet theo- 
rem (|2.7|) to show that Tot(S) = holds for surfaces obtained by a smooth 
local deformation of a plane. Thus t[il)a] can be made arbitrarily small by 
choosing a small enough. Since we want to make the form negative, we have 
to modify the trial function ip^j further in analogy with [|GJ1] . 



(d) To this aim we pick j G C^{Ao x (—a, a)) and set 9 := j'^{H — k\)iIj„. 
From (|2.14|) - (|2.15|) and the fact that the scaling acts out of the support of the 



localization function j, we immediately get the following explicit expression 

/ 2 \ 2 Xu — M 
e(g,«)=j(g,«)^vr^^J ^-^-^^^^-^sin.,n 

3 

.2 /2 



= j{q,u) 7r(-j (^(74)3 siuKiM. 

Notice that by construction the function B does not depend on a. It is 
non-zero as an element of L {VlQ^G'^cP'qdu) for a non-zero j unless G is 
independent of m. In view of (|2.1(]| ), however, the last named situation occurs 
only if K, M are zero identically on the whole surface which is impossible 
because S is not a plane by assumption. 

Since both ip^ and 6 belong to Q{H), we have 

t[^P^ + eQ]= t[^P^] + 2e\\jiH - /«?)^.f + eH[Q] . 

For all sufficiently small negative e the sum of the last two terms is negative, 
and the above arguments shows that we can choose a so that t[ilj„ + e 0] < 0; 
recall that the second term on the right side is independent of a. | 

Remark 3.3 Notice that the choice of the Macdonald function KQ{r) for 
the mollifier (f in the part (b) is not indispensable. One can modify this 
part of the proof, e.g., by using e~'~lnr. However, the choice we made is the 
most natural in a sense, because it employs the Green's function kernel at 
zero energy. 

The obtained conclusion about the bottom of the spectrum can be com- 



bined with Proposition p.l| to get the result announced at the beginning of 
Sec. O. 



Corollary 3.4 Let Q be a curved layer built over S which is a nontriv- 
ial, local, and smooth deformation of a plane, with the half-thickness strictly 
smaller than the minimum curvature radius ofL- ci. (^.ll\ ). Then H has at 
least one bound state with energy below k\. 

13 
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